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We investigate temporal behavior of probability density functions (pdfs) of paradigmatic jump- 
type and continuous processes that, under confining regimes, share common heavy-tailed asymptotic 
(target) pdfs. Namely, we have shown that under suitable confinement conditions, the ordinary 
Fokker-Planck equation may generate non-Gaussian heavy-tailed pdfs (like e.g. Cauchy or more 
general Levy stable distribution) in its long time asymptotics. For diffusion-type processes, our 
main focus is on their transient regimes and specifically the crossover features, when initially infinite 
number of the pdf moments drops down to a few or none at all. The time-dependence of the variance 
(if in existence), ~ f with < 7 < 2, in principle may be interpreted as a signature of sub-, normal 
or super-diffusive behavior under confining conditions; the exponent 7 is generically well defined in 
substantial periods of time. However, there is no indication of any universal time rate hierarchy, 
due to a proper choice of the driver and/or external potential. 
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I. INTRODUCTION 

Among a large variety of random walk and primordial 
noise models/choices, a distinguished role is played by 
heavy-tailed symmetric Levy-stable distributions, com- 
monly thought to be an exclusive support for so-called 
Levy flights. That view is founded on the concept of in- 
dependent and identically-distributed random variables 
of the jump-type and a broad class of infinitely divisi- 
ble distributions, that are admitted by generalizing the 
familiar central limit theorem, as encoded in the Levy- 
Khintchinc formula. Stable densities generically have no 
second moments. 

The response of Levy noise to external potentials is 
most often quantified by means of the Langevin equa- 
tion with an additive Levy (stable) driver. The affil- 
iated distribution functions (pdfs) obey so-called frac- 
tional Fokker-Planck equations where, in order to han- 
dle heavy-tailed pdfs, the (Brownian case) second spa- 
tial derivative is replaced by the fractional one of order 
< n < 2, c.f. (Htm and references therein. Under 
confining conditions, the second moments of resultant 
pdfs may exist and then the temporal behavior of the 
variance may be employed to quantify (on suitable time 
scales) sub-, or super-diffusive features of the underlying 
jump- type dynamics. 

In our previous papers 8 10] we have noticed that 
in addition to standard Langevin equation based meth- 
ods, an alternative modeling approach is worth investiga- 
tion. It is based on the concept of the Levy-Schrodinger 
semigroup-driven dynamics, [3 [Hj]. Contrary to the fa- 
miliar (in the context of the Brownian motion) mapping 
of the Fokker-Planck equation into the Hamiltonian dy- 
namical system, the non-Gaussian case makes a distinc- 
tion between these two dynamical patterns of behavior. 
They are inequivalent. 

Our departure point was, and still remains, the 
"stochastic targeting" (also named "reverse engineer- 
ing") strategy, [l|: given an invariant pdf p*(x), de- 
sign a stochastic jump-type process for which that pres- 



elected density is a unique asymptotic target. In case of 
Langevin-driven processes, the basic reconstruction goal 
amounts to deducing the drift function of the process. 

In Refs. [1, 0, given the very same p*(x), we have 
addressed the existence issue of a semigroup-driven dy- 
namics (e.g. the fractional version of the generalized 
diffusion equation), which relies on the existence of a 

semigroup potential V(x) = —A (| l^\^/ 2 p 1 J 2 ){x) / p 1 J 2 {x). 
Since /?* (x) is presumed to be shared with the Langevin- 
driven Levy process, we know that it has a non-Gibbsian 
functional form, [l[. Therefore standard thermalization 
and (ultimate) thermal equilibrium concepts are invalid 
in the Levy context. 

In Ref. [l(| we have relaxed the common pdf constraint 
and addressed a fully fledged reconstruction problem for 
the semigroup dynamics: given an invariant pdf, iden- 
tify the semigroup-driven Levy process for which the pre- 
scribed pdf p* may stand for a unique asymptotic one. 
Since, to this end, there is no need to invoke the Langevin 
connection, the ensuing non-Gibbsian obstacle does not 
appear anymore. 

Indeed, under new premises, asymptotic pdfs in the 
Gibbs form are admissible and a class of jump-type pro- 
cesses, non-trivially responding to environmental inho- 
mogeneities, becomes largely extended to pdfs that are 
definitely related to Gibbsian thermal equilibria, (Toj . 

The last observation suggests a possibility of a major 
recasting of the original "reverse engineering" problem 
of Ref. [1] which was designed to handle jump- type pro- 
cesses only. Namely, while before we have set common 
invariant pdfs for various jump-type processes, presently 
we shall consider invariant pdfs that are shared by, seem- 
ingly disparate, jump-type (discontinuous) and diffusion- 
type (continuous) processes. Invariant pdfs that show up 
a power-law behavior (being e.g. of the inverse polyno- 
mial form) are here allowed and thus an issue of a proper 
thermalization framework (Gibbsian equilibria), that en- 
compasses heavy-tailed distributions, reappears again. 

An issue disregarded in the past [l3[ is that generic 
Levy pdfs (like e.g. the familiar Cauchy distribution), 
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plainly against casual views, may be embedded in suit- 
able " exponential families" of pdfs fl4T4l6l | . A careful ex- 
ploitation of standard (Shannon) entropy extremum prin- 
ciples, allows to single out a concrete Levy pdf 
as a specific exponential family member, at a uniquely 
denned inverse temperature value. 

Anticipating further discussion, let us point out that 
our major observations are quite general and refer to a 
broad class of pdfs that are associated with symmetric 
stable noises and their perturbed (confining regime) ver- 
sions. It is only an analytical and numerical tractabil- 
ity reason, that makes us mostly to refer to the Cauchy 
driver in the present paper. 

We note that in the heavy-tailed asymptotic regime 
of diffusion-type processes, the process remains continu- 
ous and there is not (albeit appealing) "switch" to any 
jumping scenario. For all times we deal with a diffusion 
proceeding in conservative force fields, whose pdf asymp- 
totic features ultimately appear to mimic those normally 
attributed to jump-type processes. 

In a slightly careless manner one may think of a transi- 
tion " from a diffusive motion to the Levy flight behavior" . 
The corresponding transient phenomena (that have been 
literally interpreted as a transition from a diffusive to 
jump -typ e motion) were experimentally recorded in the 



past [171 ]. The diffusion vs jump scenario interpretation 



issue has appeared as well in the semiclassical description 
of so-called optical lattices, see e.g. [18]- 24]. 

Apart from an obvious possibility to regard the Levy 
driver (and thus to invoke the corresponding non- 
Gaussian probability distribution) as the major random 
displacement mechanism, it is perhaps less obvious that 
the Wiener driver actually may give rise to heavy-tailed 
pdfs as well. Then we encounter the dynamically gener- 
ated intermediate and/or transitional regimes, where the 
number of the pdf moments asymptotically drops down 
to a finite number or none at all. 

To elucidate these points, here we employ as a toy 
model (albeit directly related to the previously men- 
tioned optical lattice issue) a one-parameter family of 



Cauchy power pdfs, [14] : 
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a > 1/2 (1) 



as a reference exponential family of pdfs (see e.g. (14] ) 
that comprises the classic Cauchy distribution as its 
member. This familiy is an exponential one due to a 
trivial transformatin into the Gibbs-loking function with 
a logarthmic exponent. (Actually, for each "canonical" 
Levy pdf an analogous exponential family embedding can 
be accomplished.) 

The above Cauchy family arises naturally via the stan- 
dard maximum entropy principle. Namely, one seeks 
an extremum of the (dimensionless) Shannon entropy 
S(p) = — J phipdx of a continuous probability distribu- 
tion p(x) under the constraint that the expectation value 
< ln(l + 




FIG. 1: The members of Cauchy pdf family for 1 < n < 8. 



If we regard a logarithmic function a ln(l + x 2 ) as an 
external force potential, with —2ax/{\ + x 2 ) being in- 
terpreted as the forward drift in the Langevin equation, 
then we end up with a familiar Fokker-Planck equation 
(Wiener driver in action). 

A functional form of the force term is identical to that 
of the physically motivated " cooling force" in optical lat- 
tice discussions 



21] . With that restoring force, we are 



> takes an a priori prescribed value, 14] 



capable of generating all members of the above Cauchy 
family as target pdfs in the large time asymptotic of 
diffusion-type processes. This linear Fokker-Planck equa- 
tion option has appeared in the literature before, [2l|,|22j, 
also in the attempt to interpret the Cauchy family as a 
family of Tsallis distribution functions [25[ . 

All members of the Cauchy family can be achieved 
as asymptotic targets for jump-type processes. Both 
the Levy-Langevin and Levy semigroup driven dynamics 
may be employed to this end. In turn, the same "target- 
ing" is valid if we resort to diffusion-type processes. 

An important intrinsic property of the Cauchy family 
has to do with a number of moments of those pdfs. For a 
while, let us consider a parameter a in Eq. |T]) to assume 
integer values only: a — ¥ n > 1. Then one immediately 
observes that the Cauchy hierarchy of pdfs, c.f. Fig. [U 
involves a monotonically growing number (as parameter 
n grows from n = 1 to infinity) of moments of the proba- 
bility distribution, beginning from none at all for n = 1. 
The growing number of moments in existence amounts 
to an improvement of the pdf localization and entails in- 
creasing strength of confinement of Cauchy jumps (e.g. 
flights), c.f. also [1,0, if the Cauchy driver is assumed 
to be in action. 

An appealing feature of the Cauchy family is that, even 
if looking exceptional in many respects, it has received 
an ample literature coverage not only in the semiclassical 
analysis of optical lattices [lj|-[24]]. In below, we shall 
elaborate a bit on another physically important context 
of Tsallis entropies and related (via entropy maximum 
principles) Tsallis pdfs. 
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Indeed, Cauchy family pdfs were discussed in search 
for a statistical and thermal (equilibrium) foundation 
for Tsallis entropies and associated with them probabil- 
ity distributions. That includes asymptotic properties 
of various Fokker-Planck equations, nonlinear being in- 
cluded, as a dynamical justification for them, see e.g. 
Refs. [26]-{3||. From a thermal equilibrium point of view, 
in reference to heavy-tailed pdfs, our approach appears 
to provide an alternative proposal to that originally for- 
mulated in Refs. [H H, [29| and [H, HI]. 

Subsequently we shall discuss in some detail various 
random motion scenarios that may generate the above 
Cauchy family of pdfs in their large time asymptotics. 
After identifying the involved Cauchy semigroup dynam- 
ics, we shall demonstrate that diffusion-type processes as 
well generate heavy-tailed (Cauchy, to be specific) pdfs, 
in affinity with the case of jump- type processes described 
by fractional Fokker-Planck equations. 

Somewhat surprisingly, it turns out that the conven- 
tional diffusive dynamics may be regarded as an alterna- 
tive (to the fractional one) dynamical model of random 
behavior that produces non-Gaussian pdfs. This includes 
heavy tailed ones, that are casually thought to arise ex- 
clusively in relation to " free" (we stress the role of inde- 
pendent, identically distributed random displacements) 
jump-type random motions. 

In fact, we demonstrate the validity of the diffusive re- 
laxation scenario from an initial d or a Gaussian function 
to virtually any member of Cauchy power pdf hierarchy 
(fTj). In the diffusive dynamics, the resultant pdfs appear 
to have the Gibbs - Boltzmann form (as should be for 
the exponential family of pdfs). All those observations 
readily extend to a broad variety of non-Gaussian pdfs 
that have been invented in connection with the concept 
of confined Levy flights. 

Our discussion uncovers an interplay between Levy and 
diffusion-type processes that has been left basically ig- 
nored in the current " stochastic modeling" literature. We 
point out that, in principle, one may even not be able 
to discriminate between these two classes (continuous vs 
discontinuous) of random processes, if in the vicinity of 
respective (common for both) equilibrium states. Such 
observation may be particularly useful to avoid erroneous 
interpretation, in cases when an experimentally founded 
data analysis indicates a transient dynamics with a (pos- 
sibly unusual) crossover behavior, like e.g. that reported 
in Ref. [13 . 

In the present paper Langevin noises (Levy and Wiener 
driver alike) are considered as additive. For a complete- 
ness of exposition, we point out that there are reports 
in the literature that Langevin equations with the multi- 
plicative (Wiener) noise do give rise to heavy-tailed dis- 
tributions, see e.g. Ref. [37j . 



II. ENTROPY EXTREMUM PRINCIPLES 
A. Shannon vs Tsallis entropies 

An influential Ref. [l3[ appears to have harmed an 
open-minded approach to extremum entropy principles, 
with regard to heavy-tailed probability distribution func- 
tions. Its main message has been often uncritically re- 
peated: "the derivation of distributions with inverse 
power tails from a maximum entropy formalism would be 
a consequence only of an unconventional auxiliary condi- 
tion that involves the specification of the average value 
of a complicated function" . This statement has been am- 
plified at the end of Section 2 of Ref. [l3| : " It is difficult 
to imagine that anyone in an a priori manner would in- 
troduce a set of auxiliary conditions that could yield the 
logarithmic term that appears in the (Shannon) entropy 
function associated with the Levy distribution." In Ref. 
[26| . the pertinent formula (their Eq. (7)) is called the 
" ad hoc constraint" . 

A possible way out of this apparent difficulty has been 
proposed in the past. It amounts to abandoning the 
standard Shannon entropy usage in extremum principles 
and deviate towards non-extensive generalizations of the 
notion of entropy (effectively one ends up with a non- 
extensive thermostatistics), so that the heavy-tailed dis- 
tribution would "more or less naturally" arise. That has 
been suggested in Refs. [26|, l28|, |29j and followed in the 
interpretation of experimental data in Refs. (l7j|-[24l|. 

We wish to demonstrate that the non-extensive pro- 
posal is not the only admissible route, if one looks for 
mechanisms that generate heavy-tailed target pdfs. The 
conventional Shannon entropy strategy appears to work 
satisfactorily as well. 

If one seeks for an extremum of the Shannon entropy 
S(p) — — fc_B J p\xi(ap)dx of a continuous probability 
distribution (fcs is Boltzmann constant and parameter 
(j > determines the characteristic length of a sys- 
tem), under the constraint that the expectation value 
J x 2 p(x)dx = a 2 [2q |. an automatic outcome is a Gauss 
density p*(x) = ((3/tt) 1 / 2 exp(— (3x 2 ), where an (implicit) 
Lagrange multiplier f3 takes the value j3 — l/2a 2 . Finally, 
the presupposed Gibbs thermal equilibrium condition im- 
poses p = i/k B T, [mm. 

Let us indicate that for generic Gibbsian densities 
p*(x) = (l/Z) exp(—f3V(x), where V(x) stands for an ex- 
ternal potential and Z is a normalization factor, the pre- 
vious moment constraint corresponding to the Lagrange 
multiplier P may readily be generalized to the form 
(V) = const whose equivalent form is (hip*) = const'. 

So called "ad hoc" constraints of Refs. [l3[ and [26j j . 
albeit with no reference to any effective potential, nor ex- 
ponential (Gibbs-looking) recasting of the involved pdfs, 
are nothing more than the demand 

(lnp*) = J p(x)ln p*{x)dx = const (2) 

where p* stands for any a priori pre-selected Levy stable 
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pdf. C.f. Eq. (21) in Ref. [Tp and Eq. (7) in Ref. 
|26[ . We emphasize that the pertinent Levy distributions 
are "free noise" models and have nothing to do with the 
notion of confined Levy flights. 

The proposal of Refs. [H, [H, Hil amounts to consid- 
ering the new entropy function 



S q [p\ = 



1 - / d(x/a)[ap(x)] q 



with a real parameter q. An optimization [2j 
able likelihood function, under a constraint 



d{x/a)x 2 [ap{x)] q 



(3) 



of a suit- 



results in an extremal pdf: 
1 



p q {x) 



[1 - 0(1 - q)x 2 ] 



211/1-9 



where 



Z q = 



0(9-1) 



1/2 



r(i/( g -i)) 
r((3-«)/2(«-i)) 



(4) 



(5) 



(6) 



and we need 1 < q < 3 to secure convergence of the 
normalization integral. 

We note that, in view of a direct q — > 1 connection 
with the Shannon entropy (and the Gibbs density), an 
interpretation /3 = 1/fcgT is enforced. That, in turn, 
has been a starting point in Refs. [25|, [H| to deduce 
the Tsallis pdf in an asymptotic regime of a well defined 
Fokker-Planck dynamics. The resultant invariant pdf has 
the Gibbs form p q {x) — exp[— (3V(x)]/Z q , provided the 
"external force" potential reads: 



V{x) = 



1 



ln[l + P(q - l)x 2 



(7) 



Let us redefine the involved constants. Namely, after 
[2lj (we modify Eqs. (4) therein), let us set 



a 

2D' 



2D 



(8) 



That results in the V(x) 
effectively turn over into: 



1 



er 2 [l + (x/a) 2 ] and p q (x) 



P q {x) = -y [i + Way] 



2l-/3a 2 



(9) 



which differs from Eq. (fT]) by a trivial replacement of the 
exponent a by /3a 2 and an explicit usage of dimensionless 
argument x/a. Clearly, if j3 plays the role of 1/&bT, 
the signature of thermal equilibrium needs to be directly 
transferrable to the exponent a as well. 



B. Shannon entropy extremum for random systems 
in logarithmic potentials 

Cauchy family pdfs are labeled by a continuous pa- 
rameter a > 1/2, whose physical meaning seems to be 
obscure. In addition to considerations of the previous 
subsection we shall give more arguments to the contrary. 
Even without explicitly relying on the optical lattice con- 
texts, [TJ|-[I1. 

At this point we invoke a classification of maximum 
entropy principles (MEP) as given in Ref. [T3|. Let us 
look for pdfs that derive from so-called first inverse MEP: 
given a pdf p(x), choose an appropriate set of constraints 
such that p(x) is obtained if Shannon measure of entropy 
is maximized (strictly speaking, extremized) subject to 
those constraints. 

Namely if a system evolves in a potential V(x) (at the 
moment, we consider a coordinate x to be dimensionless), 
we can introduce the following functional 

/oo />oo 
V(x)p(x)dx — / p(x) \n[p(x)]dx . 
-oo J — oo 

( 10 ) 

The first term comprises the mean value of a poten- 
tial. A constant A is (as yet physically unidentified) 
Lagrange multiplier, which takes care of aforementioned 
constraints. The second term stands for Shannon en- 
tropy of a continuous (dimensionless) pdf p(x). An 
extremum of the functional L{p(x)} can be found by 
means of standard variational arguments and gives rise 
to the following general form of an extremizing pdf p*(x): 
p*{x) = Cexp(— XV (x)) which, if regarded as the Gibbs- 
Boltzmann pdf, implies that the parameter A can be in- 
terpreted as inverse temperature, A = (fc^T) -1 , at which 
a state of equilibrium (asymptotic pdf) is reached by a 
random dynamical system in a confining potential V{x). 

A deceivingly simple question has been posed in chap. 
8.2.4 of Ref. [ljj]. Having dimensionless logarithmic po- 
tential V(x) = ln(l + x 2 ), one should begin with evalu- 
ating a mean value U = (V) = V(x)p(x)dx. Next 
one needs to show that only if this particular value is 
prescribed in the above MEP procedure, Cauchy distri- 
bution will ultimately arise. Additionally, one should 
answer what kind of distribution would arise if any other 
positive expectation value is chosen. The answer proves 
not to be that straightforward and we shall analyze this 
issue below. 

To handle the problem we admit all pdfs p(x) for which 
the mean value (ln(l + x 2 )) exists, i.e. takes whatever 
finite positive value. Then, we adopt the previous vari- 
ational procedure with the use of Lagrange multipliers. 
This procedure shows that what we extremize is not the 
(Shannon) entropy itself, but a functional T with a clear 
thermodynamic connotation (Helmholtz free energy ana- 
log): 



$(x) = aV(x) +\np(x) -> T = ($> = a (V) 



S{p). 
(11) 
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Here S(p) = — (Inp) and a is a Lagrange multiplier. 
From now on, we consider an exponent a in the form 



a = e /(k B T) 



(12) 



where eo is a characteristic energy scale of a system. Note 
that here we encounter a dimensionless version of a fa- 
miliar formula F = U — TS, relating the Helmholtz free 
energy F, internal energy U and entropy S of a random 
dynamical system. 

The extremum condition 8J : {p)/8p = yields an ex- 
tremizing pdf in the form p a (x) — (1/Z a ) (l+x 2 )~ a , pro- 
vided the normalization factor Z a = (1 + x 2 )~ a dx 
exists. It turns out that the integral can be evaluated ex- 
plicitly in terms of V - functions for all a > 1/2, ending 
up at the previously introduced one-parameter Cauchy 
family of pdfs (fTJ). 

To complete an extremum procedure we can in princi- 
ple deduce a numerical value of the (Lagrange) parameter 
a, by resorting to our assumption that the mean value 
(V) a has actually been a priori fixed at a concrete value. 
This route is not at all straightforward. 

To identify the values of above a, we need an explicit 
expression for the mean value 



U a = (V) Q 



T(a) 



^T(a - 1/2) J_ x (1 



ln(l + x 2 ) 



r 2\a 



dx. (13) 



It turns out that it can be given in terms of the digamma 
function ip(x) = d(hiT) / dx: 



2tt 



sin(27ra) 



+ tl>(l — a) — ip 



a > 



1 



(14) 



This function is divergent at a = 1 /2 (see also below) and 
decays monotonously at large a. This decay is conrollcd 
by an asymptotic expansion 



1 3 1 

a ~ 2a + 8a 2 + 4a 3 



(15) 



The decay of U a at large a obeys the inverse power law. 
It follows that the expansion (fT5j) actually gives a very 
good approximation of U a for a > 3. 

We note that, apparently, Eq. (fT4")l involves another 
divergence problem, if we choose integer a. This obsta- 
cle can be circumvented by transforming Eq. (|14p to an 
equivalent form that has no (effectively removable) diver- 
gencies. Namely, we get 



U a = — 7rtan7ra + ijj(a) — if) 



(16) 



and the tangent contribution vanishes for integer a. On 
the other hand, this expression shows that the divergence 
of U a at a — > 1/2 originates from the first term in (fT6|) . as 
4> functions have finite values at this point. Near a = 1/2 
the first term of Eq. (fT6|) diverges as (a — l/2) — 1 . 

With an explicit expression for Cauchy family pdfs in 
hands, we readily evaluate Shannon entropy to obtain 



p a {x) In p a (x)dx = In Z a + a 1l a . (17) 
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FIG. 2: S a for Cauchy family (curve 1) and its asymptotic 



expansion at large a (curve 2) 
also shown (curve 3). 



Sa for Gaussian family is 



Then, the (as yet dimensionless) Helmholtz free energy 
T„ reads 



In Z n 



(18) 



with a being the dimensionless analog of the inverse tem- 
perature. 

We note, that in view of the divergence of Z a , both the 
Shannon entropy and the Helmholtz free energy (likewise 
Ua) cease to exist at a = 1/2. We plot S a as a function of 
a in Fig[2l It is seen that entropy monotonously decays 
for a > 1/2 and for larger values of a. An asymptotic 
expansion of the entropy shows logarithmic plus inverse 
power signatures 



, a \ 3 
In- + — , 

7T/ 4a 8a 2 



- ? + ... (19) 



These series are shown along with the entropy in Fig. [2] 
As a grows, the number of moments of respective pdfs 
increases. That allows to expect that an "almost Gaus- 
sian" behavior should be displayed by a > 1 members 
of Cauchy family. This is indeed the case as discussed in 
Section II.G. 



C. Thermalization via Fokker-Planck dynamics 

Let us now show that the variational principle ex- 
plicitly identifies an equilibrium solution of the Fokker- 
Planck equation for standard Smoluchowski diffusion 
processes. To address our thermalization issue correctly, 
we now use dimensional units. The Fokker - Planck equa- 
tion that drives an initial probability density p(x, t = 0) 
to its final (equilibrium) form p(x,t — > oo) reads 



d t p = DAp - V • (b(x)p). 



(20) 



Here, the drift field b{x) is time-independent and conser- 
vative, b(x) = — W(x)/(rn7) (V(x) is a potential, while 
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FIG. 3: Time evolution of pdf's p(x,t) for Smoluchowski processes in logarithmic potential ln(l + x 2 ). The initial (i = 0) pdf 
is set to be a Gaussian with height 25 and half- width ~ 1CP 3 . The first depicted stage of evolution corresponds to t = 0.01. 
Target pdfs are the members of Cauchy family for a = 1 (panel (a)), 2 (panel (b)) and 3 (panel (c)) respectively. 



m is a mass and 7 is a reciprocal relaxation time of a sys- 
tem). We keep in mind that p and bp vanish at spatial 
infinities or other integration interval borders. 

If Einstein fluctuation-dissipation relation D = 
k B T/mri holds, the equation (j20|) can be identically 
rewritten in the form dtp = V[ / oV4']/(m7), where 

# = V + k B T\np (21) 

whose mean value is indeed the Helmholtz free energy of 
random motion 

F = (*) = U - TS. (22) 

Here the (Gibbs) entropy reads S — k B S, while an inter- 
nal energy is U = (V) . In view of assumed boundary re- 
strictions at spatial infinities, we have F = —(7717) <^w 2 ) < 
, where v = —W^/(inj). Hence, F decreases as a func- 
tion of time towards its minimum F*, or remains con- 
stant. 

Let us consider the stationary (large time asymptotic) 
regime associated with an invariant density p* (c.f. Ref. 
39] for an extended discussion of that issue). Then, 
dtp = and we have VvE^p*] = Cp* (C is arbitrary con- 
stant) which yields p* = (1/Z) exp[— V/k B T\. Therefore, 
at equilibrium: 

** = V + fcsTTnp, =S> = -k B TlnZ = F*, (23) 

to be compared with Eq. (IT51) . Here, the partition func- 
tion equals Z = J exp(—V/k B T)dx, provided that the 
integral is convergent. Since Z = exp(— F*/k B T) we 
have recast p*(x) in the familiar Gibbs-Boltzmann form 
p* = exp[(F* - V)/k B T]. 

On physical grounds, V(x) carries dimensions of en- 
ergy. Therefore to establish a physically justifiable ther- 
modynamic picture of Smoluchowski diffusion processes, 
relaxing to Cauchy family pdfs in their large time asymp- 
totics, we need to assume that logarithmic potentials 
V(x) = ln(l + x 2 ) are dimensionally scaled to the form 
V{x) = 6qV where eo is an arbitrary constant with phys- 
ical dimensions of energy. 

By employing (l/eo)V(x) = ln(l + x 2 ), we can re- 
cast previous variational arguments (MEP procedure) 



in terms of dimensional thermodynamical functions. 
Namely, in view of k B T^>{x) — V(x) + k B T lnp(a;) we 
have an obvious transformation of Eqs. (fTTj) and (fl8|) into 
Eqs. (JUJ and fl23|) respectively with k B T${x) = V(x). 

The above dimensional arguments tell us that in con- 
fining logarithmic potentials, Cauchy family of pdfs can 
be regarded as a one-parameter family of equilibrium 
pdfs, where the reservoir temperature T enters through 
the exponent a. Proceeding in this vein, we note that 
eo should be regarded as a characteristic energy (energy 
scale) of the considered random system. 

We observe that a — > 00 corresponds to T — > i.e. 
a maximal localization (Dirac delta limit) of the corre- 
sponding pdf. The opposite limiting case a — > 1/2 looks 
interesting. Namely, we have S(a — > 1/2) — > 00. To 
grasp the meaning of this limiting regime, we rewrite 
a = 1/2 in the form k B T = 2eo. Accordingly, the tem- 
perature scale, within which our system may at all be set 
at thermal equilibrium, is bounded: < k B T < 2eo. For 
temperatures exceeding the upper bound T max = 2eo/k B 
no thermal equilibrium is possible in the presence of 
(weakly, i.e. weaker then, e.g., V(x) ~ x 2 ) confining 
logarithmic potentials V(x) — eo ln(l + x 2 ). The case of 
a = 1 i.e. k B T — e corresponds to Cauchy density. 

For clarity of presentation (dynamical interpolation 
scenarios between initial Gaussian and equilibrium 
Cauchy-type pdfs do not seem to have ever been ex- 
plicitly considered in the literature), in Fig. [3] we plot 
various stages of the diffusive Fokker - Planck dynamics 
for processes that all have been started from a narrow 
Gaussian. The resultant (large time asymptotic) equi- 
librium pdfs are members of the Cauchy family, labeled 
respectively by a — 1, 2 and 3. 

III. CAUCHY FAMILY TARGETING: 
TRANSIENT DYNAMICS AND RELAXATION 
PATTERNS 

To make the paper self-contained, here we describe the 
dynamics of probability distributions p(x,t), whose time 
evolution is started from the Dirac delta-type initial data 
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(actually, a very narrow Gaussian), but might be driven 
by three different stochastic mechanisms. Apart form the 
common initial data choice, we shall demonstrate that a 
common target (invariant) pdf to which the pertinent 
processes relax in the large time asymptotic, is shared 
by: 

(i) diffusion-type process with a well defined drift func- 
tion bdrift(%) and the Wiener driver in action, 

(ii) Langevin-driven jump type process, with a suitable 
drift function bj ump (x) and the Cauchy driver in action, 

(iii) Levy-Schrodinger semigroup-driven jump type pro- 
cess of Ref. [8], with an (effective) semigroup potential 
V{x) and the Cauchy driver implicit. 

The above three patterns of dynamical behavior will be 
demonstrated by assuming that asymptotic (target) pdfs 
are members of the Cauchy family 



P*n 



(l + x 2 Y 



An — 



(n-1)! 



VHXn-1/2)' 



(24) 



It turns out that for our purposes it suffices to investigate 
the above dynamics for 1 < n < 4 only, the dynamics for 
higher n being a repetition of the established patterns of 
behavior. 

Since we require processes (i)-(iii) to share an asymp- 
totic pdf p m for each choice of n, the semigroup po- 
tentials (and drift functions, if applicable) will certainly 
differ from case to case. An additional technical input 
that needs to be mentioned in view of a common initial 
pdf (actually, an " almost" Dirac delta) for all considered 
processes, is that we refer to the well developed theory of 
so-called Schrodin ger b oundary data and stochastic inter- 
polation problem, [ill. H^|. That allows to devise Marko- 
vian dynamics scenarios which interpolate between a pri- 
ori prescribed initial and terminal (eventually those refer- 
ring to the large time asymptotic and relaxation regimes) 
pdf data. 



A. Cauchy target, 

Here and subsequently, for notational and computa- 
tional simplicity, we shall get rid of all freely adjustable 
parameters (like e.g. /3, 7, see e.g. @) and set them 
equal one. The familiar Cauchy pdf: 



1 



p*i 



7r(l 



(25) 



is known to be a stationary pdf of the Ornstein- 
Uhlenbeck-Cauchy (OUC) stochastic process, (3. [T^|. 



1. Cauchy driver: jump-type process 

Following a widely accepted reasoning we take for 
granted that the Langevin equation, with additive de- 
terministic and Levy noise terms, gives rise to the frac- 
tional Fokker-Planck equation, whose form faithfully par- 
allels the Brownian version: x — b(x) + A^(t) =>• dtp = 



—V{bp) — A|AJ^/ 2 /3. All notations are consistent with 
those of Refs. [1, S HI]- To deal with the Cauchy driver, 
we set fi = l. 

Cauchy pdf (|25p is a stationary pdf of the Langevin- 
driven fractional Fokker-Planck equation, with noise in- 
tensity A = 1 

d t p = -V(bp)-\A\ 1 / 2 p. (26) 

Since the process is by construction relaxing to an invari- 
ant pdf p* (x) , the drift function is defined as follows 

1 



b(x) 



p*{x) 



(WpJ(x)dx. (27) 



Accordingly, by selecting the pdf (1231) , to which the OUC 
process actually relaxes, we identify the corresponding 
drift function as bij ump (x) = —x. 



2. Semigroup dynamics with Cauchy driver 

In Ref. we have investigated in some detail the 
Cauchy-Schrodinger semigroup-driven random motion 
scenario (e.g. semigroup dynamics with Cauchy driver). 
In this theoretical framework, the major ingredient of 
the formalism is the (pseudo-differential, self-adjoint in 
a suitable Hilbert space) Hamiltonian operator = 
A I A| M / 2 + V, that gives rise (presently, we set fi = 1 = A) 
to a fractional analog of the familiar generalized diffusion 
equation, c.f. [9j], 

|A|i/2^ 



d t ^ = -|A| i/z * - V*. 
The dynamics of the related pdf 

p(x,t) = V(x,t)pl /2 (x) 



(28) 



(29) 

is fully determined by Eq. (1281) , provided p* stands for an 
a priori given terminal (asymptotic) probability density 
of the involved stochastic process. 

The effective potential V{x) for the semigroup-driven 
jump- type process reads 

IVl 1 / 2 

V(x) = -^-. (30) 
p* 

Presently, we identify with p*\. A functional form of 
the corresponding V(x) has been found in Ref. Q: 



Vi(z) = - 



x y/ a(x) + X 



y/^(x) a(x) yfu~J~x) - x 
x 2 , c.f. also Fig. 1 in Ref. [8j 



(31) 



where a(x) =14 

The temporal behavior of the Langevin-driven and 
semigroup driven process with p*i(x) as the asymptotic 
pdf, has been visualized earlier in Figs 1 and 2 of Ref. 
jlOj. Clearly, the semigroup potential V\{x) is completely 
divorced from the harmonic potential V(x) — x 2 /2 — > 
— W = —x = bi t j ump (x), appropriate for the Langevin- 
driven OUC process. We observe that p*i(x) cannot be 
rewritten in the Gibbs form (1/Z) exp[— V(x)], in con- 
formity with arguments of Ref. [l| . 
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3. Wiener driver: diffusion-type process 

The previous two dynamical patterns of behavior seem 
to be natural in connection with the Cauchy density, 
since both refer to the Cauchy noise (e.g. driver) and 
the Cauchy operator |V| = | Al 1 / 2 as the generator of 
random motion. Therefore, it seems to be far from be- 
ing obvious that we can address the previous stochastic 
interpolation issue in terms of diffusion-type processes, 
hence with the Wiener driver in action. 

For a diffusion process that is governed by the standard 
Fokker-Planck equation 



9tP= ^Ap 



V (b ■ p) 



(32) 



with a certain {x) as its stationary solution, at least on 
formal grounds we can identify the drift of the process 
b(x) as 



1. 



°diff 



-V In 



(33) 



Consequently, upon assuming that bdiff(x) = — VV*, we 
have p*(x) = {1/Z) exp[— V*(x)//c,bT], where 1/Z is a 
normalization constant. This has the familiar Gibbs- 
Boltzmann form, although V*(x) surely has nothing in 
common with V(x) appropriate for the Cauchy-Langevin 
dynamics. 

Indeed, under the diffusion (i.e. Langevin-driven) pro- 
cess premises, 



x = b{x,t) + A(t) 



(34) 



where (A(s)) = 0, (A(s)A(s')) = S(s - s'), a demand 
that an asymptotic pdf is in the Cauchy form Eq. (f2"5"|) . 
enforces the drift function 



bi, di ff{x) 



2x 



1 



(35) 



which is a gradient field with the potential V*i(x) = 
ln(l+.x 2 ). These expressions have been substituted to the 
corresponding numerical routines to obtain a computer- 
assisted picture of the time evolution pi(x,t) and get 
confirmed that such diffusion processes are consistent. 

A comparison of the above three different (Cauchy- 
Langevin, Cauchy semigroup and Wiener-Langevin) 
stochastic interpolation scenarios, all started form iden- 
tical initial data (Dirac delta approximation by a Gaus- 
sian) data and all terminating (in the large time asymp- 
totic) at common for all Cauchy pdf, is provided in Fig. 2) 

The diffusion-type dynamics is faster than that of the 
Langevin-driven jump-type process. On the other hand, 
to achieve the target pdf, the semigroup dynamics needs 
almost twice a time (t = 15) as that (t = 8) needed for 
Cauchy - Langevin relaxation. In Fig. [4j we have indi- 
cated time instants at which the simulated pdf, within 
the resolution quality of the figures, cannot be distin- 
guished from the target pdf p*\. 






FIG. 4: Time evolution of pdf 's p(x, t) for the Cauchy- 
Langevin dynamics (panel (a)), Cauchy-semigroup-induced 
evolution (panel (b)) and the Wiener-Langevin process (panel 
(c)). The common target pdf is the Cauchy density, while the 
initial t = pdf is set to be a Gaussian with height 25 and 
half- width ~ 1CP 3 . The first depicted stage of evolution corre- 
sponds to t = 0.01. The time rate hierarchy seems to be set: 
diffusion being fastest, next Levy-Langevin and semigroup- 
driven evolutions being slower than previous two. However 
the outcome is not universal, as our subsequent discussion 
will show. 



The above relaxation pattern is generic for Cauchy 
family of pdfs. However, the resulting "speed" hierarchy 
is not generic for stochastic processes at hand. Below 
we shall demonstrate that time rates of relaxation for 
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our three exemplary processes (i.e. jump-type Langevin, 
jump-type semigroup, diffusion-type) do not follow any 
definite hierarchy. 



B. Second order target, p„2- 

With the basic notation specified, we can proceed in 
quick steps. Let us consider the n — 2 pdf of the Cauchy 
hierarchy 



P*2 



r ; 1 + x 2 ) 2 



(36) 



If regarded as an invariant density of the Cauchy- 
Langevin stochastic process, this pdf gives rise to the 
following drift function 



^2, jump 



(x) = -^-3) 



(37) 



The semigroup-driven process, relaxing to p*2(x) in the 
large time asymptotics, is determined by specifying the 
semigroup potential, [|| c.f. Fig. 1 there-in: 



V 2 (x) 



x 1 - 1 

X 2 + 1' 



(38) 



The affiliated diffusion process (Wiener driver), with the 
very same asymptotic pdf has the drift function 



b 2 .d l ff{x) = 2b hdiff (x) 



Ax 



1 + x 1 



(39) 



At this point, we note a recursive formula b n>( nft(x) = 
nbi^ di ff{x) and V m (x) = n ln(l + x 2 ). 

Qualitatively, the temporal behavior of the three in- 
volved stochastic processes mimics the pattern of n — 1, 
according to Fig. HI except that with the growth of n the 
speed with which the processes relax increases. 



where a(x) — 1+x 2 . Like the previous two semigroup po- 
tentials, that were bounded from below and above, the 
present one, although unbounded from above, still ful- 
fills minimal requirements set in the general theoretical 
framework of Ref. 12]. It is bounded from below by 
^3(0) = — 4/-7T and for |x| — > 00 we have Vs(x) w 2\x\/ir. 

The diffusion process with (x) as its asymptotic tar- 
get is characterized by the drift function &3,di//(a;) = 
-6x/(l + x 2 ). 



D. Fourth order target, p*4. 

Although we know that nothing illuminating might 
happen in the qualitative picture of the time evolution, 
if compared to the previous cases, for completeness we 
provide drift functions and the semigroup potential, that 
can be derived if n=4 pdf is to be asymptotic target of, 
respectively: Cauchy-Langevin, Cauchy semigroup and 
diffusive motions: 
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5tt(1 + x 2 ) 4 ' 
The Cauchy-Langevin drift presently reads 



(43) 



b^ump(x) = -y^ 35 + 35x 2 + 21x 4 + 5x 6 ) , (44) 

the Cauchy-Schrodinger semigroup potential appears in 
the form 



Va{x) 



x 4 + 6x 2 - 3 
2(x 2 + l) ' 



(45) 



while the drift of an affiliated diffusion process is 
b4,diff(x) = -8x/(l + x 2 ). 



E. Tail (large \x\) regularities 



C. Third order target, p*3. 
The third, n=3 member, of the pdfs family reads 

p * 3 = (40) 

It is an invariant probability density of the Cauchy- 
Langevin process with the drift function 



b3,jump(x) = - — (15 + lOx +3x ). 



(41) 



An effective potential for the Cauchy-Schrodinger dy- 
namics can be evaluated by means of the Cauchy prin- 
cipal value integrals. Following methods of Ref. [8j we 
have arrived at its explicit functional form 



V 3 (sO 



2(x 2 -2) 3x , Ja(x) + x 

y ' + —— In v ; v ' 

yfa{x) a{x) ^Jaixj-x 



(42) 



It is easily seen from the above equations that as \x\ — > 
00, we encounter some regularities. Namely, pdfs in the 
Cauchy hierarchy for large |a;| behave like 



P*n(x) « „ 



(46) 



Cauchy-Langevin drifts have large |x| asymptotics of the 
form 

bn,jump(x -> OO) S3 X 2 "" 1 . (47) 

Cauchy-Schrodinger semigroup potentials behave like 

21n|x| 1 



Vi(N->c») 



x\ \x\ 

— I |0 



V 2 (\X\ ~>00) S3 1 = |x| U , 
V 3 (|JC| OO) S3 \x\, 

V 4 (\x\ -»■ oo) S3 \x\ 2 , 



(48) 

(49) 
(50) 
(51) 
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so that for arbitrary n 

V n (\x\ -> oo) « \x\ n - 2 . (52) 

The forward drifts of a diffusion - type processes, 

for which Cauchy family pdfs (|24l) are their respective 
asymptotic targets, read 

2n 

K,diff( x ) ~ — -■ ( 53 ) 



F. A cross-over between exponential and power 
law behavior 



of this quantity in the course of the OUC process, with 
that of the second moment for a diffusion-type probabil- 
ity distribution p(x,t) having p*i(x) ([2"S")) as the target 
(asymptotic) pdf, see Fig. [5]i. 



G. Interpolation between Cauchy and Gaussian 
distributions 

Since, within the Cauchy hierarchy, the number of mo- 
ments increases with n, we can justifiably ask for its link 
with the Gaussian pdf. At least, for large n. 

Let us rewrite Eq. ([TJ) in the form 



We have mentioned before that the Cauchy hierarchy 
is ordered with respect to an overall number of moments 
in existence, of the probability distribution. Namely, for 
n = 1 even the first moment does not exist. In principle 
(if Cauchy principal value of the corresponding integral 
is considered) we can admit the existence of the first mo- 
ment. But still, the second moment does not exist. 

For n = 2, we have first three moments, the third one 
exists conditionally if the above Cauchy principal value 
is considered. Ultimately, each p m has exactly 2n — 1 
moments for each n > 1. 

On the other hand, in our numerical procedures, the 
nonlinear (in a sense that corresponding drift function is 
not linear) diffusion - type motion is started from a very 
narrow Gaussian. As indicated in the caption to Fig. 21 
its height (maximum) is ~ 25 while half-width ~ 10~ 3 . 
This pdf certainly admits all moments. 

In the course of the diffusion - type evolution with a 
target pdf belonging to the Cauchy family, quite apart 
from incompatible microscopic mechanisms (jumps vs 
continuous paths) there must arise a definite cross-over 
from Gaussian to Cauchy hierarchy behavior. The lat- 
ter hierarchy exhibits the power law asymptotics to be 
set against Gaussian (stronger then exponential) decay. 
Therefore higher moments of the initial pdf, while evolv- 
ing diffusively, must consecutively disappear in the course 
of time, if any Cauchy density is to be a target pdf. 

To visualize this intriguing crossover between Gaussian 
and Cauchy family behaviors, we have numerically com- 
pared the time development of largest finite moment for 
a given n (which is (2n — 2)-th moment of Cauchy family 
(|2"4")l pdfs) with that of a smallest nonexistent moment 
(which is 2n-th moment of the same family). This evo- 
lution has been generated by Cauchy - Langevin driven 
diffusion-type p(x, t) beginning (at t = 0) at above very 
narrow Gaussian and ending (at t — > oo) at p m (|2"4")l . 
This has been done for n = 2 — 4 and is reported in 
Fig. El 

The special case n = 1 has been considered separately. 
Namely following our earlier results, Ref. [lfj, we con- 
sider the half-width at half-maximum of function (1251) in- 
stead of the above largest finite moment. Then, similar 
to the above n > 1 case, we compare the time evolution 



p a (x) 



(1 + x 2 ) 



- = exp(\np a ) , 



(54) 



where lnp Q = lnj4 a — aln(l + x 2 ). 

We observe that the limiting procedure a — > oo en- 
forces strong localization properties of the pdf so that 
only a small vicinity of x = matters. Indeed, at large 
a, the pdf tails become very steep and do not play a deci- 
sive role. This motivates the following series expansions 



ln(l + y) 

A a (a — > oo) 



Accordingly: 



p a (x,a > 1) 



y 2 , y 3 

v 

y 2 3 

(i 3 

7T 8\/7ra 



V - + ..., y = x 2 , (55) 
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128a Jna 



.(56) 



7 



7r 128av / 7ra 



x 

(57) 



exp -ai + ex- 



it is obvious that leading terms (in the large a regime) 
coincide with a perfect Gaussian 



p a (x,a > 1) 



— exp(— ax ). 

7T 



(58) 



We can easily check that function (|58jl is normalized. 
Moreover, Eq. (|57p sets a "bridge" between the Gaussian 
and large a members of the Cauchy family. The account 
of a sufficiently large number of terms in the series (|57|) . 
permits to control an accuracy of the " Gaussian approxi- 
mation" of Cauchy family pdfs with large a. Ultimately, 
since a ~ > oo stands for a sequential approximation of 
the Dirac delta functional, our limiting procedure within 
the Cauchy family provides another sequential approxi- 
mation of the Dirac delta. 
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(a) Half- width and < X 2 > (divergent) for p„i (b) < X 2 > and < X 4 > (divergent) for p*2 




(c) < X 4 > and < X 6 > (divergent) for p„z (d) < X 6 > and < X 8 > (divergent) for p»4 



FIG. 5: The time evolution of the last convergent and first divergent moment for pdfs p tl - p„4, shown on panels (a) - (d) 
respectively. Crossover time instants are shown on the corresponding panels. 



IV. GENERIC RELAXATION PATTERNS 
A. Regularity conditions in stochastic modeling 

1. Jeopardies 

The present subsection is motivated by various regular- 
ity conditions that guarantee the existence and unique- 
ness of solutions, both for stochastic differential equa- 
tions and Fokker-Planck equations (fractional being in- 
cluded), see Refs. and Those mathematical 
restrictions do not seem to worry physics-oriented prac- 
titioners and we shall indicate why in some cases this 
disregard might be justified. 

A transparent picture of jeopardies to be met, comes 
within the traditional Brownian motion framework. 
There, drifts and so-called diffusion functions (in our 
case, they reduce to diffusion coefficients) need to obey 
the Lipschitz and growth bounds, [1^, |4(|. It is easy 
to generate Langevin and Fokker-Planck equations with 
drifts that violate the growth condition. Then, some care 



is necessary, because the general theory tells us that the 
stochastic process in question may explode in a finite 
time, [40| , 

Polynomial drift functions typically violate the growth 
condition. We shall analyze a possible significance of 
this violation, by employing: (i) numericallly generated 
solution of the Fokker-Planck equation, (ii) trajectory - 
wise Monte Carlo (MC) simulations of pdfs, in terms of 
Langevin equation-generated bunch (c.a. 10 5 ) of Brown- 
ian paths. 

We demonstrate that, despite the violation of the 
growth condition, no hint of a possible explosion is de- 
tected and the dynamical pattern of behavior remains 
regular up to arbitrarily large times, available in our nu- 
merical simulations, both in FPE solution and in direct 
MC simulations. Note, that in preceding discussion we 
have already invoked the FPE (both standard and frac- 
tional) with polynomial drifts. A number of examples 
can also be found in Refs. d, 0, [H|- 
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FIG. 6: Time evolution of a diffusion process with drift b(x) = —4x 3 . Panel (a) shows the pdf dynamics, obtained from 
FP equation. Panel (b) reports the comparison of path-wise MC simulation for t = 5 with analytical result for asymptotic 
pdf. Simulated pdf has been intentionally left "noisy" (i.e. created with smaller number of sample trajectories) to reveal the 
coincidence with analytical curve. Panel (c) compares the MC result at time t — 10 with the analytical result for p*2(x) = 
(2/7r)/(l + x 2 ) 2 . The simulated pdf has been intentionally left "noisy". 



2. Cubic drifts 

Let us consider a cubic function b(x) = —Ax 3 as an 
exemplary drift in the F-P equation. This function surely 
does not obey the growth conditio n |39l .l40j: \b(x) \ + 1 < 
C(l + |x|). Standard arguments pj,[39l|allow us to identify 
a stationary pdf associated with this drift function. It 
has an explicit Gibbs form = (1/Z) exp(— Vi), where 
V* = x 4 . 

To get convinced that it is actually an asymptotic pdf 
of a well defined diffusion-type process, we address the 
time evolution of p(x,t) that is started from "almost" 
Dirac delta pdf and propagated in accordance with the 
corresponding FPE. The resultant relaxation to the non- 
Gaussian equilibrium pdf p*(x) ~ exp(— x 4 ), is depicted 
in Fig. IH^i. 

This relaxation pattern has been directly path-wise 
confirmed by executing the MC simulation of c.a. 10 5 
Brownian paths, up to the time instant t = 50. This 
allowed us to check the shape and stability of the cor- 
responding pdfs. Fig. [BJ} reports a direct comparison 
between the simulated (up to time t = 5 as for t > 5 the 
corresponding pdfs are exactly the same) and analytic 
asymptotics p*(x) = 0.55163 exp(— x A ). 

Our simulation shows, that a perfect agreement of 
two different generation methods for a given pdf extends 
to (arbitrarily) large times. Each of the transitional 
pdfs shown in Fig. [SJi can be consistently reproduced 
trajectory-wise, i.e. by means the MC simulation whose 
record is stored and analyzed at each required time in- 
stant. In Fig. we have intentionally allowed the simu- 
lated asymptotic pdf to be somewhat " noisy" . That out- 
come is a consequence of not too large number of sample 
trajectories employed for the pdf approximation. Thus, 
on computer-assisted grounds, the time evolution of pdfs 
generated by FPE's with polynomial drifts, appears to 
be consistent. 



3. Cauchy family targeting 

In relation to the discussion of Section II, our simula- 
tions have shown that a proper choice of the drift function 
in the FPE, as t — > oo takes the initial (^-function into 
any Levy - stable distribution target, or the like (e.g. 
members of the Cauchy family). Accordingly, conven- 
tional diffusion-type processes may in principle give rise 
to "heavy-tailed" asymptotic pdfs. Thus, heavy tails are 
not an exclusive property of Levy processes. 

The only (and very substantial) difference between 
"normal" and fractional FPE asymptotic outcome is that 
in the fractional case Levy flights may not relax to Gibbs- 
Boltzmann pdfs, [l|. In turn, diffusion modeling with a 
suitable external forcing admits a Gibbs form of asymp- 
totic "heavy-tailed" pdfs. 

Fig. |BJ; reports an analogous (F-P dynamics vs MC 
simulation) comparative procedure, in case when both 
Cauchy and Wiener drivers are responding to different 
external forces, while giving rise to a common asymptotic 
pdf. Cauchy dynamics does not produce the Gibbs pdf. 
However, the same pdf has the Gibbs form in the diffusive 
relaxation process. 

We have directly checked that a target function be- 
longing to Cauchy family (l24l) . p^{x), is an asymptotic 
outcome of the diffusion process via path- wise simulation, 
with b2,diff{x) of Eq. (f3l?|) . We have confirmed that the 
MC procedure (with check-out times t = 0,01, t = 0, 1, 
t = 1 and t = 8) correctly reproduces the qualitative 
picture of the time evolution depicted in Fig. 0^,. The 
ultimate outcome in the vicinity of the asymptotic pdf is 
reported in Fig. [6}:. 

The very same reasoning may be adopted for exem- 
plary jump-type processes, with Cauchy driver and poly- 
nomial drift functions. We note that relaxation prop- 
erties and confinement of Levy flights in various exter- 
nal potentials has received attention in the literature, see 
e.g. references in [J-ll, For example, the Cauchy- 

Langevin dynamics has been analyzed for dtp — — |V|p — 
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FIG. 7: The relaxation to the bimodal pdf (|59|) for three different type of processes: Langevin driven (panels (a) - (c)), semigroup 
driven (panels (d) - (f)) and Wiener driven (conventional FPE diffusion) (panels (g) - (i)). The initial Dirac delta- like pdf is 
located, respectively, at x = — 1 (panels (a), (d), (g)), x — —0,5 (panels (b), (e), (h)) and x = (panels (c), (f), (i)). All three 
processes have different drift functions, corresponding to terminal pdf (|59[1 as indicated in the text. On panels (a) and (b), the 
log scale is utilized for better visualization of time evolution of initial unimodal pdf to the terminal bimodal one. 




FIG. 8: The time evolution of variances for pdfs in Fig. [7] Panels (a), (b) and (c) correspond to Langevin driven, semigroup 
driven and conventional diffusion (Wiener driven) process respectively 
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V(6/o), with bj ump (x) — — a; 3 . The corresponding non- 
Gibbsian invariant pdf reads p*(x) = 1/tt(x 4: — x 2 + 1). 
0,1- 

At this point we can again invoke p*2 (%) , this time re- 
garded as an asymptotic target of the Cauchy-Langevin 
jump- type process with the drift function &2,j«mp(^) = 
— x(x 2 — 3)/8 of section II. B. Both the resultant pdf evo- 
lution and the corresponding Monte Carlo simulation, in 
terms of jump sample paths, give consistent outcomes. 

Concerning mathematical subtleties of the Cauchy 
semigroup dynamics, their wcll-dcfinitcncss comes from 
the properties of semigroup potentials. This issue has 
received due attention in Ref . [l2| , see also @, H[ . 



B. Relaxation time rates 

In recent years there has been continually growing in- 
terest in various theoretical random walk models, clearly 
motivated by the inefficiency of the standard Gaussian 
modeling paradigm. That refers to various sub-fields of 
physics, extending to chemistry, biology, biophysics and 
financial mathematics. 

The classical concept of Brownian motion pervades the 
whole theory of stochastic processes. The pdf of a homo- 
geneous Brownian motion solves a Fokker-Planck (actu- 
ally, heat) equation and remains sweeping unless confined 
in a finite region or subject to external forces. Typically 
the time-evolving pdf, that is initially concentrated at (or 
about) a point, with the flow of time takes the Gaussian 
form, whose width grows in time as t 1 / 2 . This diffusion 
processes is called the normal diffusion. 

In a broad field of anomalous diffusions (and more gen- 
eral non-Gaussian processes), whose microscopic mod- 
eling may involve both Markovian and non-Markovian 
random dynamics, another property has been discovered: 
< X 2 {t) >^ t 1 , with < 7 < 2. This dynamics is sweep- 
ing as well, hence precluding the existence of any well 
defined relaxation pattern (except for a possible relax- 
ation to a uniform distribution, if the motion is spatially 
confined). 

For a special class of non-Gaussian jump-type pro- 
cesses, first and second moments of the pdf may not exist 
at all. That happens in case Levy - stable distributions. 
Due to long tails of the pdf, in the least their second (and 
all higher) moments are nonexistent. The existence of the 
first moment is granted only for a suitable subclass. 

In the presence of external potentials the confinement 
(taming) of Levy flights may occur, and quite in affin- 
ity with standard diffusion processes in external force 
fields, asymptotic invariant pdfs may be approached. 
The speed (time rate) of the corresponding relaxation 
processes is worth addressing and set against that aris- 
ing in the sweeping motion. 

Somewhat intriguing point in the anomalous transport 
is whether, often used terms like " subdiffusion" or "su- 
perdiffusion" may at all have meaning in connection with 
Levy flights, under confining (i.e. relaxation) regimes. 



Clearly, there always appears a transitional period dur- 
ing which the diffusion, jump-type or semigroup-driven 
process may be analyzed in terms of a time dependent 
variance (unless non-existent). This property is regarded 
as generic in relaxing to equilibrium diffusion processes. 

Our purpose is to check whether one may expect any 
regularity in the time evolution of variances for diffusion- 
type and jump-type (Langevin and semigroup-driven) 
processes. That could possibly give useful hints towards 
their hierarchical classification with respect to relaxation 
time rates. May we talk about "slow", "fast" or "super- 
fast" processes at all ? 

Such hypothesis might look plausible on the basis 
of our previous analysis of relaxation scenarios whose 
asymptotic targets were pdfs from the Cauchy family, 
Eq. ([1]). However, the situation is not that simple and 
appears to be incongruent with naive expectations. This 
issue we shall discuss below. 



1. Relaxation to the Cauchy pdf p*i 

Our first observation pertains to the Ornstein- 
Uhlenbeck-Cauchy process (OUC, Langevin-driven) and 
its semigroup-driven relative with a common target pdf 
p*i, [_}• ^ i s seen form Fig. [H that both the transitional 
behavior and the time needed to reach (almost, within 
the figure resolution limits) the asymptotic pdf, indicate 
that the semigroup dynamics is running somewhat slower 
as compared to Cauchy-Langevin-driven dynamics. In- 
deed, the transitional period of motion admits a direct 
analysis of the half- widths (second moments of p*\ are 
non existent, see above) time rate in terms of respective 
exponents of t 1 . Namely, for the proper OUC process, 
the transitional exponent equals 7 = 0.58, while in the 
semigroup-driven case, we have 7 = 0.45. 

Both processes are of the jump- type and have no sec- 
ond moments. We note that the Cauchy-Langevin driv- 
ing sets the p(x, t) in the vicinity of the invariant one af- 
ter time t = 8, while the semigroup-driving needs t = 15 
to the same end. The diffusion scenario, according to 
Fig. QJ: sets at p*i after time t = 5. 

2. Relaxation to the bimodal pdf 

Now we discuss various aspects of relaxation time rates 
for diffusion-type, Cauchy-Langevin driven and Cauchy 
semigroup-driven processes. All of them are started from 
"almost" Dirac delta pdfs (albeit localized about three 
different initial points, and all terminated at an asymp- 
totic bimodal pdf, [l(| 

_ 2a{a 2 + b 2 ) 1 

p^x) _ - ________ 

(59) 

Here, parameters a and b are, respectively, real and 
imaginary parts of complex roots of the cubic equa- 
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tion z 3 + z— 1/4 = 0. Their approximate values are 
a ~ 0.118366 and b ~ 1.0208. 

This exemplary bimodal pdf is an asymptotic distribu- 
tion of the Cauchy-Langevin driven dynamics of p(x, t) 
with the drift bj ump — — W*(x), where V*(x) = (x 2 — l) 2 . 
The Cauchy-Langevin evolution of p(x,t), that is started 
from three different locations x — 0, x — —0, 5 and 
x = — 1 in the double well potential V*, is depicted in 
Fig. panels (a)-(c). 

The data collected for Fig(7](a)-(c) allow to deduce the 
time-dependence of variances, in the transitional regime. 
They display ~ t 1 behavior which depends on the initial 
pdf location data. Namely, for x — — 1 and x — —0, 5 we 
have respective exponents 7 = 0.53, while x = yields 
7 = 1. 

The Cauchy semigroup-induced evolution of p(x,t), 
has been simulated under the very same (as in the 
Cauchy-Langevin case, Fig 13 (a)-(c)) initial data. The 
corresponding patterns of behavior are depicted in FigJTl 
panels (d)-(f). In Fig. [S]we report the time evolution of 
respective variances. The ~ t 1 behavior in the transi- 
tional regime is characterized by exponents: 7 = 1 for 
x = —1, 7 = 0.53 for x = —0.5 and 7 = 0.2 for x = 0. 

For completeness we have performed analogous sim- 
ulations with the Wiener driver in action, i.e. for the 
conventional diffusion-type process that interpolates be- 
tween the common (for all three types of processes) ini- 
tial data and the terminal bimodal pdf, see e.g. also 
Ref. [43j]. The results are shown on panels (g) - (i) of 
FigEJ In Fig. [3] we also report the temporal behavior of 
related variances. The 7 exponents in the transitional 
regime read: for x = —1.5 and x = —1 we get 7 = 0.75, 
x = —0.5 corresponds to 7 = 0.5, while x = to 7 = 1.2. 

The preceding analysis shows that standard classifica- 
tion of anomalous diffusion processes as " subdiffusions" 
or " superdiffusions" on the basis of the exponent 7, is 
invalid in the presence of external forces. Moreover, this 
exponent strongly depends on both the particular local- 
ization of the initial pdf (effectively, an initial position of 
a fictitious particle), with respect to the corresponding 
potential profile, and on the potential curvature. 



V. CONCLUSIONS 

The main message of the present paper is that, un- 
der suitable confinement conditions, the ordinary Fokker- 
Planck equation may generate a non-Gaussian heavy- 
tailed pdf (like e.g. Cauchy or more general Levy stable 
distribution) as its long time asymptotics. That implies 
a continuous interpolation between an initial highly lo- 
calized pdf (like Dirac 6 - function) having all moments 
and the terminal heavy - tailed pdf, with only few (or 
none at all) moments in existence. 

Since it is the fractional generalization of FPE which 
is customarily invoked to generate the same heavy-tailed 
pdfs, albeit with a very different choice of drift functions, 
we have uncovered an unexpected link between non- 



Gaussian jump-type processes (inherent in the fractional 
FPE derivation) and ordinary diffusion processes that are 
based on the Gaussian (Wiener driver) paradigm. 

In the present paper we have paid attention not only 
to random processes that have found their place in the 
literature, all stemming from different forms of Langevin 
modeling (here, only additive noise has been considered), 
but also to alternative jump - type scenarios based on the 
concept of the semigroup-driven dynamics, [1, [l|. This 
pattern of dynamical behavior has not been satisfactorily 
explored within the area of anomalous random transport. 

To make our findings transparent, we have undertaken 
a direct numerical verification of them, by invoking two 
kinds of modeling. The first one amounts to solving 
numerically both ordinary (with the Laplace operator) 
and fractional FPE's with the drift functions selected to 
generate, the same for both, non-Gaussian (specifically 
Cauchy family (j2"4")l ) pdfs. The second one presumes that 
one has the above drift functions at disposal. Then, we 
undertake a direct numerical modeling of random trajec- 
tories (sample paths) of the underlying stochastic pro- 
cess, both Brownian and of the jump-type. These tra- 
jectories have next been used to generate the associated 
pdfs which turn out to coincide with those obtained from 
numerical solution of both ordinary and fractional FPE's. 

We note that a traditional interpretation of random 
data relies on the existence of at least two lower mo- 
ments of involved probability distributions. Therefore, 
the case of Levy flights could have been placed within 
this interpretational paradigm only under confining con- 
ditions. In such case, in view of the existence of sec- 
ond moments, one may address their time-dependence in 
transitional regimes, when the process is yet far from its 
equilibrium pdf. The above numerical solutions of the 
FPE and of its fractional analog permit us to investigate 
the time dependence of different moments (see Fig. 8) of 
the corresponding pdfs. 

At this point we have passed to an approximation of 
the variance a 2 by the power law a 2 ~ t 7 , with the hope 
that its validity would shed some light on the problem 
of assigning terms like "subdiffusion" or "superdiffusion" 
to anomalous transport processes. Namely, the time evo- 
lution of variances is often considered as a signature of 
subdiffusive, normal or superdiffusive behavior, typically 
considered with no bounds on the duration time of the 
process (unless set be experiment). Taken literally, all 
such motions would belong to the sweeping category, 
where the variance grows indefinitely 

In this connection, our analysis shows that (i) near ini- 
tial t — > the exponent 7 depends strongly on time, (ii) 
deeply in the asymptotic t — > 00 regime 7 can be equated 
to 0, (iii) in a transitional regime (i.e. for pdfs lying some- 
where between initial and terminal ones) the exponent 7 
happens to take almost constant positive values (observe 
e.g. almost straight lines in time dependence of variances 
in the middle parts of the corresponding plots in Fig. 8). 

We have focused on possible signatures of (ab)normal 
relaxation patterns, expecting that in a transitional 
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regime we might possibly identify an undoubtful subd- 
iffusive or superdiffusive dynamics. Our answer is neg- 
ative. We have confirmed that widely used terms "sub- 
" and/or " superdiffusion" are rather vague and become 
meaningless for confined stochastic processes, likewise in 
the presence of Levy and Wiener drivers. The pertinent 
exponent value is not indicative for the stochastic process 
in question, to justify its naming (sub-, normal, super-), 
because it depends strongly on the location (position) of 
the initial pdf relative to the drift potential profile, and 
on its curvature. 

It appears that the above "anomaly" in the descrip- 
tion of random motions appears to have much deeper 
foundations. One can alternatively use incompatible mi- 
croscopic mechanisms to model an engineered asymptotic 
approach to the same a priori prescribed non-Gaussian 
target pdf. 

If those pdfs are an outcome of a statistical analysis 
of experimental data, clearly a distinction between dif- 
fusive and jumping patterns of dynamical behavior ap- 
pears not to be that sharp as commonly expected. On the 



other hand, a possible discrimination tool, e.g. the single- 
particle/molecule experiments in the nano- or mesoscopic 
domain should be set under scrutiny, since microscopic 
features of motion may not be verifiable at current lev- 
els of precision/resolution in collecting the data. The 
specific stochastic model that would seemingly fit to the 
data, might possibly be not more than a lucky guess, to 
be invalidated in a more thorough analysis. 

Another message worth spelling is that appar- 
ently divorced dynamical mechanisms (Langevin-Cauchy, 
Cauchy semigroup and diffusion-type process) may share 
common asymptotic pdfs. Contrary to the wide-spread 
belief, heavy-tailed probability distributions, casually as- 
sociated with jump-type processes, may as well have 
Gaussian (e.g. diffusive) origins. That is possible if drift 
functions refer to properly functioning potentials, which 
we attribute to a cumulative effect of inhomogeneities 
of the environment). Their role would be to attenuate 
the casually considered counterbalancing of the Wiener 
driver, "normally" necessary to eliminate heavy tails of 
the resultant pdf. 
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